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Abstract 

One-loop matching of heavy-light currents is carried out for a highly improved lattice action, 
including the effects of mixings with dimension 4 0{1/M) and 0{a) operators. We use the NRQCD 
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improved glue action. As part of the matching procedure we also present results for the NRQCD 
self energy and for massless Asqtad quark wavefunction renormalization with improved glue. 
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I. INTRODUCTION 



The importance of approaching the chiral hmit in a controlled way in lattice simulations 
of heavy-light systems has become increasingly evident in recent years. Progress in reducing 
statistical and discretization errors has resulted in chiral extrapolation uncertainties, to- 
gether with operator matching errors, dominating the final total error in heavy-light meson 
decay constant and form factor calculations. This has motivated the HPQCD collabora- 
tion to initiate studies of heavy- light systems using improved staggered/naive light quarks, 
taking advantage of the good chiral properties of such actions which allow one to go down 
to much smaller quark masses than has been possible in the past [1, 2, 3, 4]. Heavy-hght 
simulations are now being carried out with light quark masses as low as m strange I'^i ^-nd 
contact has been established with chiral perturbation theory predictions. 

An important ingredient in all heavy meson decay constant and form factor calculations is 
the matching of the heavy-light currents used in the simulations to their continuum QCD 
counterparts. The highly improved heavy-light actions introduced above necessitate a new 
round of matching calculations. In this article we report on the one-loop perturbative 
matching of the temporal component of the heavy-light axial current, ^40, with NRQCD 
heavy quarks and massless Asqtad naive quarks. These matching coefficients are directly 
relevant for our ongoing Jb, /bs ^^'^ fos calculations on the MILC dynamical configurations 
and have already been applied in the results quoted in [1, 2, 3]. Due to the chiral symmetry 
of naive quarks, matching coefficients for the vector current are identical to those of the 
axial current; thus the results presented here can also be applied to form factor calculations 
[5]. Our matching calculations include contributions from 1/M current corrections and an 
0{a as) discretization correction. The matched heavy- light current is then correct through 
0{as), 0{a tts), 0{as/{aM)) and 0{as I^qcd/ M). Further corrections would come in at 
0{al), 0{KIcd/M^) and ^(a^ a,). 

In the next section we list the quark and ghie actions employed and give a brief discussion of 
our calculational methods. Two independent strategies were adopted and the results tested 
against each other. The wavefunction renormalization constant Zq for massless Asqtad 
fermions, which is used in later sections on matching coefficients, is presented in section 3. 
Section 4 describes one-loop self-energy corrections for NRQCD heavy quarks. Both the 
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heavy quark wavefunction renormalization Zq and the mass renormahzation Zu enter into 
the current matchings. Section 5 presents the full mixing and matching calculations for the 
NRQCD/Asqtad heavy-light currents. We give tables of results for a range of heavy quark 
masses. We add several appendices with calculational details, such as a list of Feynman 
rules. 



II. THE LATTICE ACTIONS AND CALCULATIONAL STRATEGIES 



The Glue Action 



The tree-level tadpole and 0{a}) improved glue action is given by [6], 



So = -ti E 11 

a;, ti>v y 



(1) 



ut 12 12 j ' 
with 

P^. = -^Re (Tr{U^{x)U,{x + a^)Ul{x + aMi^)}) , (2) 

P/.^. = ^Re (Tt{U^{x)U^{x + a^)U^{x + 2a^)Ul{x + + a^)Ul{x + aMi^)]) ■ (3) 

/3 = 2Nc/g^ and uq is the tadpole improvement factor [7]. We carry out perturbative 
calculations in general covariant gauge and add a gauge fixing term, 

n 2 



^5,/ - ^ E 



(4) 



with d^j^A^j^{x) = A^{x + — An{x — a^/2) and ^ the gauge parameter. We work in 

both the Feynman and the Landau gauges and verify that gauge invariant quantities are 
independent of ^. In one-loop matching calculations of current operators, the only ingredient 
required from the gauge action is the tree-level gauge propagator. The form of the free gauge 
propagator for improved gauge actions has been discussed in several articles [6, 8, 9]. We 
follow reference [9] and invert the momentum space free improved gauge action once and for 
all using Mathematica. One ends up with a closed form analytic expression for the gauge 
propagator which is a 4 x 4 matrix in Lorentz space. The interested reader is referred to 
Appendix A of [9] for more details. 
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B. The Light Quark Action 



The "Asqtad" action is the most successful quark action to date for simulating light 
dynamical and valence quarks [10, 11, 12]. The MILC and HPQCD collaborations are using 
the MILC dynamical configurations to study light-hght, heavy-light and heavy- heavy physics 
[13, 14, 15, 16, 17] [2, 3, 4]. This action has been designed to drastically reduce the taste- 
symmetry breaking effects of conventional staggered fermions and also has all other 0{a^) 
discretization errors removed. As discussed in references [11] [1] there are two equivalent 
ways to write down the Asqtad action, either in terms of four component "naive" fermions or 
in terms of one component "staggered" fields. We prefer to use the naive fermion approach. 
We will use Feynman rules for four component fields and our heavy-light currents will be 
constructed out of improved naive light quark fields. 

Three new ingredients are necessary to go from the unimproved staggered/naive action to 
the Asqtad action, fattening of links, the so-called "Lepage correction" and the Naik term. 
One starts from the unimproved naive action. 



with 



^(x) , (5) 



V^^{x) = [U^{x)^{x + a^) -Ulix- a^)^{x - a^)] , (6) 

and Hermitian Euclidean 7-matrices obeying {7^,7;^} = 25^,^. The link fattening and the 
Lepage terms are incorporated by replacing the link variable Uf^ in the covariant derivative 
by an updated variable (the corresponding covariant derivative will be denoted V^). 

y;(a;) ^ V,{x) - Y: ^^U,{x) (7) 
V,{x) ^ n 1 + ^ U,ix) (8) 

pj^H \ / symmetrized 

Vfj,{x) is the fattened (taste violation suppressing) link and the second term in eq.(7) is the 
Lepage term that removes a low momentum O(a^) error. The derivatives, Vf, and Vf,'*^^\ 
that act on link matrix variables are defined as {/i ^ u), 

-VlU^ix) = -L\u^{x)U,{x + a,)Ul{x + a,) 
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-V^'^'^ U,{x) = \ \U^{x) U,{x + a^) Ulix + a,) 

lO Uq 

--U,{x). 



(9) 



(10) 



The Naik correction adds a three-link term to the modified covariant derivative so that the 
final form of the Asqtad action takes on the form 

1 /_, 1. 

with 



^Asqtad 



X] 



" 6 



m 



(11) 



^3— link 



8 1 





3 



tadpole improved 

UUU ^{x + 3a^) - U^U^U^ ^{x - 3a^) 



U "i/ix + a^) - U^ *(x - 



(12) 



Feynman rules for the action (11) have been written down in concise form by Q. Mason [18]. 
We verified the rules for the one-gluon emission vertices and for the subset of two-gluon 
emission vertices relevant to the present calculations, namely the subset symmetric in the 
two gluons. These Feynman rules are summarized in the appendix. 



C. The Heavy Quark Action 

We use the NRQCD action improved through (9(1/M^) and O(a^), and which also in- 
cludes the leading relativistic 0{1/M^) correction [19]. This action is currently being used in 
simulations of heavy-heavy and heavy-light systems on the MILC dynamical configurations. 
It has been discussed in many previous publications, and hence we will be brief here. In 
terms of the two-component Pauli spinor one has. 




5 



Ho is the nonrelativistic kinetic energy operator, 



and SH includes relativistic and finite-lattice-spacing corrections, 

aSH = -ci ^ S + C2 2 ( V • E - E • V ) 

-C3 — ^—2 <T • (V X E - E X V) 

8(aMo)' 



1 ~ A(4) (A(2))2 

— C4 — cr • B -I- C5 — ^ — Cr n . (15) 

^2(aMo) '24(aMo) ' 16n(aMo)' ^ ^ 



All derivatives are tadpole improved and. 
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A(2) = ^vf, A(4) = ^vf (16) 

Vfe = Vfc - ^Vl^^ . (17) 

Precise definitions of the V^*'*, i — 2,3, 4, and of the improved E and B field operators are 
given, for instance, most recently in Appendix B of [1] (we note, however, a factor of 2 error 
in eq.(B4) of this reference, which is corrected above in eq.(9)). Feynman rules for simpler 
versions of NRQCD actions have appeared in [20, 21, 22]. We have generalized them to 
rules for the more highly improved action (13) - (15) studied here. The new Feynman rules 
are given in an appendix. 



D. Calculational Strategies 

We have employed two independent approaches to carrying out the one-loop matching 
calculations with the above actions and used them as checks against each other. In one 
method we make extensive use of Mathematica to multiply propagators and vertices, to carry 
out the Dirac algebra, and, where necessary, to take derivatives with respect to external 
momenta. The resulting expressions are then converted into FORTRAN and fed into a 
FORTRAN VEGAS code for the numerical integration over internal momenta. 

In many instances the Mathematica expressions can be kept quite simple and general, since 
details of Feynman rules are only needed at the VEGAS stage. For instance, in calculations 
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that do not involve any derivatives with respect to external momenta, the only information 
the Mathematica code requires is the formal structure of the vertices in the Asqtad and 
NRQCD actions. For the emission of a gluon with polarization //, the Asqtad action has the 
general vertex, 

V{A,)^Y.^,,,^^ (18) 

and the NRQCD action leads to 

Vh{A,_,) = wh^fi + ^wh^j aj (19) 
j 

aj = Pauh matrices (see end of Appendix A.l for a 4 component version of (19)). The 
dependence of the "^"'8 on the incoming and outgoing fcrmion momenta is specified only 
in the VEGAS code through subroutines that code up the Feynman rules. 

Mathematica is also very useful when one needs to take derivatives with respect to external 
momenta. One inserts the Feynman rule expression for the Wfj_^i, or whu^u of interest and lets 
Mathematica take derivatives before setting external momenta to zero. 

Our second method is based on a C++ code developed by C. Morningstar for matching cal- 
culations with NRQCD/clover quark actions and which was generously made available to us. 
We modified the original code to replace clover by the Asqtad hght quark action and also to 
accommodate improved glue with non-diagonal (in Lorentz space) gauge propagators. The 
(7+"*" code handles derivatives via "automatic differentiation." A C"*""*" class is defined which 
carries information not just about a function, but also about its first couple of derivatives 
through a Taylor series expansion. When two such class instances are multiplied, for ex- 
ample, derivatives of the product are calculated automatically via the chain rule as part of 
the definition of an overloaded multiplication operation. We believe our two methods are 
sufficiently independent of each other so that they provide good checks of our results. 

III. Zq FOR MASSLESS ASQTAD QUARKS 

Perturbative calculations with the one-component (improved staggered) version of the 
Asqtad action and with unimproved Wilson glue have been performed recently in [23] for 
the matching of light-light bilinear and four-fermion operators. Calculations using the four- 
component (improved naive fermion) Asqtad action and improved glue have been carried out 
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for bilinear and four-fermion operators [24] and for mass renormalization [25] . For the match- 
ing of heavy-hght currents the hght quark wavef unction renormahzation is required. This 
quantity has been derived aheady with improved glue using "twisted" boundary conditions 

[26]. We have repeated the calculation with a gluon mass IR regulator, the IR regulator 
we use throughout this article. is relevant for many other light-light and heavy-light 
perturbative calculations. We list our results here for the case of massless light quarks. 

= 1 + a, {C\'' + C,] + 0{al) (20) 

= ± [1 + (e - 1)] ln(a2A2) (21) 

OTT 

where A is the gluon mass. The gluon mass dependence cancels, of course, upon taking 
the difference between continuum and lattice one-loop coefficients. The IR finite Cg has 
contributions from the regular rainbow diagram, the tadpole diagram, and from the 
tadpole improvement procedure. We write 

= + + Cf (22) 

and hst results in Table I for both Feynman and Landau gauges. The uq correction contri- 
bution is given by 



4 2 



«f=-j«?'. (23) 



(2) 

with Uq defined through, 

uo = l-as 4'^ + 0{al) . (24) 
The three contributions in (23) come from the fat-link, Lepage, and Naik terms respectively. 

(2) 

The numerical value of Uq depends on how one defines uo and on the gauge action. In the 

(2) (2) 

case of improved glue one has Uq ' — 0.767 for the plaquette definition of uq and Uq ' — 0.750 
for the Landau link definition. The MILC dynamical configurations were created using 
the plaquette uq. It then makes sense to use plaquette uq in nonperturbative valence light 
quark propagators and in perturbative light quark self energy calculations. For the NRQCD 
heavy quark parts of heavy-light simulations we have adopted the Landau link uq and this 
is reflected in the heavy quark parts of the perturbative matching given below. There is 
no contradiction in using different definitions of uq in different parts of a calculation. The 
important thing is that the perturbation theory match the choices made in the numerical 
simulations. We give explicit formulas for the tadpole improvement contributions (such as 
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FIG. 1: Self-energy diagrams which generate Cq (or Cq if a hght quark hne is substituted). 

(23) and several others below) so readers should be able to adjust the numbers in our Tables 
to their uq choices. 

IV. HEAVY QUARK SELF ENERGY 

One-loop self-energy calculations have appeared already in many articles for a variety of 
NRQCD actions on the lattice. Refs. [20, 22] present results for simpler NRQCD actions 
with unimproved glue, Ref. [27] deals with highly improved NRQCD actions and Wilson 
glue, and Ref. [28] works with Symanzik improved glue but at low order in NRQCD. Here 

we summarize results for the NRQCD action of eqns. (13) - (15) and improved glue. We 
are interested in wavefunction renormalization Zq, mass renormalization Zm, and (for com- 
pleteness) the energy shift Eq. 

Zq = 1 + asiClJ' + Cq] + 0{al) (25) 
Zm ^ I + asCu + 0{a1) (26) 
aEo ^ asCEo + 0{al) (27) 



where. 



^l-[-2 + {i- 1)] Ha'X') . (28) 



Just as for light quarks, the one-loop heavy quark self energy Q;sS(p) gets contributions 
from the two diagrams in Fig. 1 and from tadpole improvement. If one defines 



fio = S(0) (29) 

(30) 



0,1 = -i 



^2 = 2{aM) 



9S 



(31) 

Pn=0 
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where stands for the heavy quark external momentum measured in units of the inverse 
lattice spacing (S is also taken to be dimensionless) and M is the heavy quark pole mass, 
then, 



Ceo = ~ ^0 

Cq = fio + 

Cm — ~ ^1 ■ 



IR finite 



(32) 
(33) 
(34) 



Results for the one-loop coefficients Cxx arc tabulated in Table 11 for several values of (oMq) 
and n. To the order we are working there is no need to distinguish between the pole mass 
M and the bare mass Mq in one-loop coefficients. Since simulations are carried out at fixed 
(aMo) it is much more convenient to express everything in terms of the bare mass. The 
results in Table II were obtained using the Landau link uq, with the explicit formulas given 

by, 



— 



-1 - 



1 + 



C5 



(aMo) 
3 



2 (aMo) 

C5 



+ 



Cl 



(aMo; 



+ 



C6 



3 2n(aMo)2^ 

Cl , Cq 



,.(2) 



+ 



(2) 



(35) 
(36) 



(aMo)3 2n(aMo)2^ 

There are no tadpole improvement corrections to Cq since such corrections cancel between 
and ill, as do the contributions from the tadpole diagram. Wherever possible we give 
analytic formulas as functions of the coefficients q in the NRQCD action eq.(15), however 
numerical results such as in Table II are always given with Ci — 1. 

Zm and Eq are gauge invariant, however Zq is both gauge variant and IR divergent. We 
discuss in an appendix how we handle IR divergent terms such as in (25) and isolate the IR 
finite contribution. 



V. MATCHING OF THE HEAVY-LIGHT CURRENT 

The basic formalism for perturbative matching of heavy-light currents has been developed 
already for NRQCD/clover currents [21] and can be taken over without any modification for 
the case of improved naive light quarks. Again we will be brief and refer the reader to the 
earher articles for details. The following three currents are needed in the NRQCD effective 
theory to match the temporal component of the vector or axial vector currents to full QCD 
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through 0{as/M). 



4'\x) = q{x)roQ{x), 



-1 



4'\x) = 



2(aMo 
-1 

2 (aMu 



-q{x) ro7-vg(x), 



(37) 
(38) 

(39) 



The Q fields are four component Dirac spinors with upper two components given by the 
NRQCD (p field and vanishing lower components, and Fq = 70 or 7570. The V in Jq^ and 
Jq'^ is the same as in (6). The relation between matrix element of Aq in full QCD to those 
of the effective theory currents can be written as, 

(0)\ 



(40) 



with 



J, 



{i),sub 







as do 4'^ (41) 

for i — 1,2. We prefer to express things in terms of the more physical matrix elements 
(J^^^'""'') which have had C»(a,/(aMo)) power law contributions subtracted out [29]. The 
coefficients Pq, pi and p2 are then given by. 



Po 



lie, 



a 



Q) 



00 



Bi — -^{Cq + Cq) — Cm 



Pi 

P2 — B2 — C02 — C12 



Coi — Cii 



(42) 



with 



Bo ^ - 

TT 
1 



Bi 
B2 



TT 

4 

TT 



In(aMo) - - 
19 

In(aMo) - - 



(43) 



Cq, Cq and Cm are the one- loop coefficients of self-energy corrections discussed in the two 
previous sections. Cm enters into pi because we have written the 1/M currents in (38) and 
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FIG. 2: One loop vertex corrections which generate the dj. 

(39) using Mo rather than the pole mass. Originally these currents are defined in terms of 
the pole mass, the mass that is common to both full QCD and the effective lattice theory. 
One then, for convenience, replaces M by ZmMq and notes that to one-loop order one can set 
Zm — 1 everywhere except in the tree-level contribution from Jq'^ , where Zm = 1 -I- q;^ Cm 
must be used. 

The in (41) and (42) come from the mixing between the three currents. They are 
calculated at one-loop from diagrams shown in Fig. 2 where one puts Jq^ at the vertex 
and projects out the tree-level expression {Jo^)tree- For instance, in the case of the Cio that 
appear in the power law subtraction (41), one is taking the matrix element of Jq^ and asking 
how much of it projects back onto Jq'^ Since we match at zero external momentum one 
has Cio — C20- As explained in references [21, 30] ^2, or more specifically C02, includes a term 
that removes an 0{a ag) discretization error from Jq'\ The matching procedure is such that 
0{as/M) and 0{aas) corrections are (and must be) made at the same time. 

We mention that some care is required in calculating , = 0, 1. There is another dimension 
4 current operator involving a time derivative acting on the light quark field q{x) which is 
equivalent to Jq via equations of motion. 

When calculating C,i2 one must also include contributions from {JQ^)tree- 

Our results for po, Pi, p2 and ^10 are given in Table III. Only pi has a Uq correction, coming 
from the contribution to from Fig. 3: 

Cit = . (45) 

We use the Landau hnk uq in calculating pi. Table III can be used in several ways. Should 
one choose not to include matrix elements of the 1/M currents Jq^^ and J^^ in one's simula- 
tion, then matching should be done with just the first term in eq.(40). In other words, one 
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FIG. 3: uq correction to pi: C|\ . 



should not include the — ag Cio {Jq ) term which we have absorbed into ( Jq ). This way 
one avoids unnecessarily introducing a lattice artifact ©(^^/(aMo)) power law term that 
can only be canceled by the matrix element {Jq'^)- Using just the first term in (40) leaves 
us with 0{Aqcd/M) and 0{aas) errors. The next step in improvement would be to use 
(^o) = (1 + CKsPo) {Jq'^) + (Jq^^'*"^)- Again there are no 0{as/{aMo)) power law contri- 
butions and errors from the heavy-light current operators come in now at 0{as Aqcd/M) 
and 0{aas). The latter two errors are removed by going to the full expression in (40). 



VI. SUMMARY 



We have performed a one- loop matching of the lattice NRQCD/Asqtad heavy-light cur- 
rent (temporal component) to its continuum QCD counterpart, correct through 0{as/M) 
for a range in heavy quark mass. One-loop results for the NRQCD heavy quark self energy 

corrections and the massless Asqtad quark wavefunction renormalization are also presented. 
We find that all the perturbative coefficients are well behaved and none of them are partic- 
ularly large. 

The matching coefficients are important ingredients in lattice investigations of B and D 
meson leptonic and semileptonic decays. The 0{a1) errors that remain after incorporating 
the one- loop results of this article will, at some point, become the dominant error in decay 
constant and form factor determinations from the lattice (this is already happening for /^^ 
and calculations [2]). Consequently, there is a need to push on to two- loop matching 
calculations. Members of the HPQCD collaboration are already engaged in higher order 
perturbative calculations with improved lattice actions [31, 32, 33]. It should be possible to 
extend those calculations to higher order perturbative matching of heavy-light currents as 
well. 
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APPENDIX A: FEYNMAN RULES 



In deriving Feynman rules for fermionic actions we rely heavily on methods developed by 
Colin Morningstar in reference [27]. The basic quantity of interest is the "^-function". 



= iV{k', k ; (gi, i/i, 61), (g2, ^2, ^2), ■ ■ ■ , (gr, J^r, W) ) 



(Al) 



which gives the vertex due to operator A for a fermion to emit r gluons of polarization 
i/i, i/2, . . . , i^r, color index bi,b2, ■ ■ ■ ,br and with momenta qi,q2, ■ ■ ■ ,qr- k' is the momentum 
of the outgoing fermion and k the momentum of the incoming fermion, k — k' + qi- 
Given such vertices for two operators A and B, one can derive the emission vertices for 
the product operator AB using the convolution theorem of Fourier transforms. The general 
expression is given in eq.(40) of [27]. Here we reproduce only the special cases of zero- one- 
and two-gluon emission vertices, which are the only ones needed in the one-loop calculations 
of this article. 



C%{k',k;{q,u,b)) 



ife^k', k; {qi, i/i, 61), (^2, 1^2, ^2)) 



Sk'M^f{k,k) ^f{k,k) 



A2 



{k',k; iq,u,b)) Cf{k,k) + 



ik',k') {q,i^,b)) 



(A2) 
(A3) 



{k',k; (gi,z/i,6i), (g2,f^2,&2)) ^f{k,k) + 
{k',k') ^f{k',k; {qi,UiM),(Q2,i^2,b2)) + 



^f(k', k' + qi; {qi, Ui, bi)) ^f{k -q2,k; {q2, 1^2, ^2)) 

(A4) 

These equations are extremely useful and also very general. The operators A or B can 
be simple operators, such as the link variable U/^ or V^^. They can also stand for more 
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complicated combinations such as F^^^ all the way to terms hke (l — or (l — Ul{t — 

1) (l-^)" ^ that appear in the NRQCD action. 

Reference [27] lists the for r = 0, 1, 2, 3, for a wide range of elementary operators, C/^, 
V^, A(^\ F^j^i, and for improved versions of these operators. Using the convolution equations 
(A2) - (A4), one can then build up the emission vertices for composite operators such as 
V • E or (A^^))^ or a string of U matrices. This is done, for instance, in the C+"^ code 
(method 2 of section IID). We have also used these equations to obtain Feynman rules once 
and for all by hand. They are summarized in the following two subsections. 



1. Feynman Rules for the NRQCD Action 

In this subsection we give the Feynman rules for the NRQCD action (13) - (15). We 
will do so in the two component language which is appropriate for self-energy and heavy- 
heavy current matching calculations. At the end of this subsection we mention the trivial 
generalization to a four component language which is more useful for heavy-light current 
matchings. 

In order to streamline the notation and also to introduce spin information, we re-express 
the relevant ^^J^ as follows. 



(r) Vl,...Vr 



ej\k',k;{qi,Ui,h),iq2,U2,b2), (g.,i^.,M)^ ^ ' {k' , k) T'K...T'^ (A5) 

The gluon momenta qi,....,qr are left implicit in this notation. The labels s = 0,1,2,3 
refer to vertices with spin structure (J2, o'x,o'y,o'z) respectively. We also rewrite the NRQCD 
action in terms of the following operators: 

0.= (i-^) (A6) 
O, = fl-^)" (A7) 



and schematically one has. 



2n 

O2U2 = O2 U} O2 (A8) 



'-NRQCD — 

0t(l - 0,02U20i)(f>. (A9) 
Feynman rules are obtained by determining, 

[0^02U20,]^:^''''■■■'''^{k',k) (AlO) 
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where we are using the notation of (A5). 
Tree- level: 



One has 



[0,]f{k\k) = 5s,oh',kF,{k) 



with 



F(k) = 



Cl 



C6 



(aMo)3 2n{aMof 



and 



[02t/2]f {k', k) = 5,,o e-* F(fc)2" 
The free heavy quark propagator is given by, 

G%\k) = {l - [0,0,u20,tlo{k,k)Y" 
= {l - e-''=°F(A;)'"Fi(A;)2}"' . 



(All) 
(A12) 

(A13) 
(A14) 

(A15) 



(A16) 



One-gluon emission: 

Using the notation of this appendix, the coefficients wh^j^^g of eq.(19) are defined (up to a 
color T-matrix which we omit) as, 

= [Oi02c/20i]il^^ {k\ k) , wh,,^ a, = [Oi02c/20i]il^; {k' , k) (A17) 

Repeated use of (A3) leads to, 

[OilSo k') [0,u2]^^' {k', k) [OilSo {k, k) (A18) 

= [Oi]f^^(A:',A:) [e-< F{k'f'' Fi{k') + e''^' F{kf'' Fi{k)) + 
F,{k')F,{k) [02U2t^^ {k\k) . (A19) 

Going further down the chain one has, 

[02[/2]i'^''=° {k\ k) = 8s,o F{k'f F{kf (-1 e-'('=o+'=o)/2 j ^^20) 
[02U2]^s^^=' {k' ,k) = [02]^^' {k' ,k) {e-<F(kT + e''^' F{kY) (A21) 
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with {k^ = ^) 



-1 



[O^ir^^ {k', k) = 5s,o ^^^^ sin(A;;) F,(A;', A;, n) 



F2{k\k,n) = < 



1 n = 1 

E?=o F{k'Y-^-^ F{ky n > 1 



and 



[Oi]'rik':k) = -^j:m]'^^^{k',k). 



(A22) 



(A23) 



(A24) 



The \bH,^l^^ {k\k) can be built up from the Fourier transforms hsted in [27] using again 
(A3). One finds, 



Cl 



+ 



^ {i^ ,k) 

C6 



2(aMo)3 4n(aMo)2 



sin(/c^) 



7=1 ^ ;=i ^ 



(A25) 



- I W^t^' ik', k) = 5s 



a 
2 



16(aMo)2 ^ 



E {sin(2A;;"~) cos(ko) [sin(k'i) - sin(ki)] rjio] 



[(^i/sli^^^ (k', k) = (5,,o -^Q^^l^^Y {sin(2/co ) cos(/c^. ) [sin{k'^) - sin(%)] 77^0} 



(A26) 



(A27) 



^ [^^3]i=- /j) = '^i E ^iJ-' cos(A;o ) [ss(/cp + ss(/c^)J sin(2/c; ) r/zo 

(A28) 

C3 



^ ['^^^3]S = iq^^Im^y cos(/Cj. ) [ss{k'i) + ss(/cO] sin(2/co ) r/jc 



(A29) 



I [<^^4]i2(A;',A;) = 4(^]^ ^ ^^j' cos(A;^- ) sin(2A;, )r]ji 



(A30) 
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C5 



6(aMo) 



sm{k^) + - sm{2k^) cos(A; • 



(A31) 



where 7]^,, — 1 for unimproved E and B fields and 77^,^ = |[5 — cos(2/c~) — cos(2/c~)] for 
improved fields. The function ss{kj) — sin{kj) for unimproved Vj and ss{kj) — |sin(A;j) — 
|sin(2A;j) for the improved Vj. 

We have coded up the one-gluon emission vertices, [02(72] g^'*'^ (/s', A;), [Oi]f^'^ {k' , k), etc. as 
functions of arbitrary incoming and outgoing momenta, k and k'. These routines are called 
repeatedly by the VEGAS code when calculating the integrand of the one-loop diagrams. 
Similar expressions are also fed into Mathematica routines, whenever one needs to take 
derivatives with respect to external momenta. 



Two-gluon emission: 

We will restrict the discussion to two-gluon emission relevant for the tadpole diagram, for 
which only a special combination of momenta are required : {k',k,qi,q2) — > {k,k,q,—q). 
One is also interested only in the spin singlet contribution. Applying (A4) to this case one 
has, 

[Oi02U20i]fl!^^'{k,k) = [Oi]fl^'''{k,k)2e-''°F{kf''Fi{k) + [O2U2& {k,k) Fi{kf + 

Fi{k) [02U2]'-^l'o{k,k + q) [Oi]2o(^ + ?,fc) + 

[O£l'o (k, k + q) [02U2]flo {k + q, k) F,{k) + 



{[0^]^^^'{k,k + q) [0,]^^>{k + q,k)) 



-i(ko+qo) 



s=0 



F{k + q) 



2n 



(A32) 



The two new ingredients are [Oil^^o'"^ (k, k, g, —q) and [02U2\flo''^ {k, k, q, —q). 
[0,]fl^,''{k,k,q,-q) = l,,^ + 5,, 5., 



+ 



2(aMo)3 4n(aMo)2 

3 U. 1 

1' ' 2 

iC2 



SijCOs{kj) 5]sin2(^) ^ isin(/ci + |) sm{kj + |-) 



16(aMo)^ 



(5^j5^,o + S^,oS^.,j) cos(% + ^) sin(g^) cos(^) rjjo 



- S^j6i,j 2 cos{kj + sin(go) cos(^) rjjo 
+ i2(aMo) ^^'^^'"'^ ^^^^^^'^ ~ cos(2%)cos2(|-) 



(A33) 
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where rjjo is the same as the r]^„ defined after eq.(A31) with 2A;^ — > q^. 

[02U2&''{k,k,q,-q) = 2e-^'^Fikr[02&{k,k,q,-q) - 5,,,5,/-e-^'^^F{kf- 

- iF{kre-^^''°^'^^ [S,,o [02]^llo{k + q,k) + S.,o [02^1^ {k, k + q) 
+ e-('=°+«°) [O2]2o^ {k, k + q) [OslSo' + 9, k) (A34) 

where 



[02]fl'o''{k,k,q,-q) = 6,,A 



^,3 



1 \^ • 

' sin(/cj + ^) sm{kj + ^) F^{k, k + q,n) 



2n{aMo) J ' 2 ' ' ^ 2 



(ASS) 



with 



,0 n = 1 

F3(A;,A; + g,n)= <^ ^ . (A36) 

Er=~i ^ + g)""'"! n > 1 

4 component formulas : 

So far we have written down Feynman rules corresponding to the NRQCD action (13) which 
is given in terms of 2 component Pauli spinors. One can equivalently write rules appropriate 
for the Q fields of (37) - (39). One needs to make the replacement 

{l2,a) ^ (/4,E) (A37) 

with Ej = diag((7i, (jj) denoting 4x4 matrices. Furthermore, the above gluon emission 
vertices should be multiplied by ^-y^. 

2. Feynman Rules for the Asqtad Action 

As mentioned already in the main text, Feynman rules for the Asqtad action for one- 
gluon emission and the most general two-gluon emission vertices have been written down by 
Q. Mason [18]. In this appendix we hst just those rules used in our calculations. 

Tree-level : 

The free Asqtad propagator is given by, 

= E7^sin(M [l + ^sin2(A;^)] + (am) | . (A38) 
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One-gluon emission : 

Carrying out the derivatives and symmetrizations of eq.(8) one ends up with a fattened hnk 
consisting of [1-hnk, 3-hnk, nonplanar 5-hnk, nonplanar 7-hnk] staples with weights [1/8, 
1/16, 1/64, 1/384] respectively [12]. The Lepage term leads to a 1-link contribution to 
with weight 3/8 and a planar 5-link contribution with weight -1/16. Adding to this the Naik 
term one ends up with. 



[Asqtad] 



[fat link] ^^^'^ + [Lepage] ^'^'^ + [Naik] 



(i)m 



i(i)/' 



(A39) 



with (we again define = 



[fat linkj^^)" {k', k) = cos(A;+) cos^ik') cos^ik;) cos^k;) 



+ Il7;^cos(A;+)sin(A;^)sin(A;^) \^ \ {cos^{k, 



[Lepage]«'^(^',A;) = -i | 7^cos(A;+) ^ sin2(2^;) 



COS' 



'iK)) 



(A40) 



- 7v cos(/c+) sin(/c^ ) sm{k^ ) cos^(/c^ ) 



(A41) 



[Naik](^^^(A;',A;) = -i <[ 7/.cos(/t+)- [cos\k+) - cos\k-) (l - 4sin2(A;+ 



(A42) 



In the above formulas, the Lorentz indices n, z/, p, a are all taken to be different from each 
other. One can read off from these equations the coefficients W/^^i, of eq.(18). 



Two-gluon emission : 

We will again present only those two-gluon emission vertices specific to the tadpole diagram 
where {k', k, qi, ^2) {k, k, q, —q). One has 



[Asqtad] 



{7.sin(fc,) [1 + 1 ([l-4cos^(fc,)][l-4cos2(|)]2 + 3 
+ E7.sin(Msin^(|)2sin2(|) 



(A43) 



[Asqtad] 



- <{7^sin(A;^) 



-^sin(|)sin(|) (2 + cos^(|) + cos^(|) 



+ 2cos^(|)cos^(|)) + ^sin(|)cos(|)sin(g.) 
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APPENDIX B: INFRARED SUBTRACTIONS 



At intermediate stages of the matching procedure several of the lattice one- loop integrals 
are IR divergent. We use a gluon mass A to regulate those integrals and extract the IR finite 
contributions in the following way (similar approaches can be found, for instance, in [34] [9]). 
If J-hik-i A, Mo) is the relevant lattice integrand, one can write. 



[Mk,X,Mo) = / [J^L{k,X,Mo) - Tsui,{k,A,X,Mo)] + R{A,X,Mo), (Bl) 

Jk Jk 

with 

R{A, A, Mo) = / J^subiK A, A, Mo) (B2) 

Jk 

and 

To be consistent with the previous Appendix, the momentum k is always measured in 
units of the inverse lattice spacing. A, A and Mq however, have dimensions of energy. 
J-'suh{k, A, A, Mq) is constructed to have the same IR behavior as the lattice integrand of 
interest, so that the first integral on the RHS of (Bl) is IR finite. A is a (ultraviolet) cutoff 
imposed on the subtraction term so that J^sub{k, A, X, Mq) — for k'^ > a^A^. We take 
A < ^ and vary A to check that the dependence on this cutoff cancels between the two 
terms in (Bl). Another criterion for choosing a suitable Tgub is that the integral in (B2) be 
easy to do and that one be able to extract the IR divergent piece, i?^^, analytically: 

i?(A, A, Mo) = i?^*"^*"(A, Mo) + i?^^(A, A, Mo) . (B4) 

Otherwise there is a lot of freedom in choosing Tsub- We give one example for the logarithmic 
IR divergence in the NRQCD wavefunction renormalization Zq 

The IR divergence in Zq comes from Vti defined in (30). A convenient choice for J- sub in 
this case is the corresponding quantity in continuum NRQCD, which must have the same IR 
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behavior as the lattice theory. One can easily convince oneself that the full continuum S is 
actually not required. The IR divergence resides in the contribution to S from a temporal 
gluon in the one-loop rainbow diagram of Fig. 1, the only contribution that survives into 
the static limit. Setting /j, = u = in the gluon propagator one finds, 



J^^,ik, A, A, Mo) = < 



(B5) 



with, 

w = -T—:rT ■ (B6) 
2 (aMo) ^ ^ 

Using 4D spherical coordinates the angular integrals in (B2) can be carried out analytically 
(we use MAPLE for this) and one is left with a ID radial integral. 



i?(A,A,Mo) 



1 r'^' dy 



471 Sn Jo {y + {aXy 



2vy5 + 1 



3/2 



+ 



2(e-l) 



y6 + 3^/^+l- + 1 



3/2 



yb 2VP + 1 



3/2 



where. 



1 



(B7) 



fB8) 



4 (aMo)2 

The IR divergent contribution to R{A, A, Mq) can be isolated straightforwardly and one has. 



j^finite 



v 



a2A2 



dy 



{y + (aA)2) 



2Vy5 + 1 



3/2 



- 2 



+ 



(e-1) 

a2A2 



y6 + 3Vy5+l - [2Vy5 + 1 



3/2 



dy 



Stt Jo {y + (aA)2) 
^l2-«-l)l4 



1/6 [2^/yb + 1 
[2 - (e - 1)] 



3/2 



+ 1 



(B9) 



aA->0 



(BIO) 



The A dependence is canceled by the first term on the RHS of (Bl). The gluon mass 
dependence is the same as in continuum QCD. 
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TABLE I: Coefficient of tlie one-loop contribution to Zq for massless Asqtad quarks. The plaquette 
definition of uq is used to implement tadpole improvement. ^ is the gauge parameter. Where no 
errors are indicated, they are of order one or less in the last digit. 



i 




Qtad 






1 


-0.798(3) 


1.600 


-1.726 


-0.924(3) 





0.000(3) 


1.310 


-1.726 


-0.416(3) 



TABLE IL Coefficient of one-loop contributions to the heavy quark self-energy. The Landau link 
definition of is used to implement tadpole improvement. ^ is the gauge parameter. Where no 
errors are indicated, they are of order one or less in the last digit. 



aMo 


n 


<^Q(e = l) 


CQ(e = 0) 


Cm 




5.40 


1 


-0.016(3) 


0.493(3) 


-0.026(4) 


0.917 


4.00 


2 


-0.142(3) 


0.368(3) 


0.082(4) 


0.850 


2.80 


2 


-0.338(3) 


0.171(3) 


0.235(4) 


0.767 


2.10 


4 


-0.539(3) 


-0.030(3) 


0.387(4) 


0.666 


1.95 


2 


-0.611(3) 


-0.105(3) 


0.421(4) 


0.689 


1.95 


4 


-0.603(3) 


-0.094(3) 


0.428(4) 


0.646 


1.60 


4 


-0.797(3) 


-0.288(3) 


0.543(4) 


0.609 


1.20 


6 


-1.184(3) 


-0.678(3) 


0.732(4) 


0.609 


1.00 


6 


-1.545(3) 


-1.037(3) 


0.859(4) 


0.758 
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TABLE IIL Matching Coefficients for the temporal components of the heavy-light axial and vector 
currents. The Landau link uq was used to evaluate pi. Where no errors are indicated, they are of 
order one or less in the last digit. 



aMo 


n 


Po 


Pi 


P2 


ClO 


5.40 


1 


0.240(2) 


0.509(4) 


-0.155(6) 


-0.095 


4.00 


2 


0.149(2) 


0.445(4) 


-0.134(6) 


-0.123 


2.80 


2 


0.043(2) 


0.382(4) 


-0.072(6) 


-0.166 


2.10 


4 


-0.044(2) 


0.346(4) 


0.002(6) 


-0.210 


1.95 


2 


-0.058(3) 


0.343(4) 


0.008(6) 


-0.222 


1.95 


4 


-0.065(3) 


0.342(4) 


0.018(6) 


-0.219 


1.60 


4 


-0.117(2) 


0.343(4) 


0.061(6) 


-0.259 


1.20 


6 


-0.175(2) 


0.385(4) 


0.107(6) 


-0.329 


1.00 


6 


-0.186(2) 


0.463(4) 


0.119(6) 


-0.378 



25 



